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PART  CXVI 
HYDRAULIC  DESIGN 

CHAPTER  6 
SURGES  IN  CANALS 


^4)1.  SCOPE 

This  chapter  is  confined  to  the  problem  of  surges  which  occur  in  navigation  canals  as  a  result 
of  lock  operations,  since  a  full  discussion  of  other  unsteady  flow  phenomena  which  may  occur  in 
open  channels  would  be  feasible  only  in  a  large  treatise. 

а.  Surge  Problems.  In  canal  design,  the  study  of  surges  is  necessary  because  of  the  following 
three  problems  which  are  given  in  their  probable  order  of  importance:  (1)  to  an  extent  governed 
largely  by  traffic  density;  (2)  surges  which  depress  the  water  level  reduce  the  effective  channel 
depth,  those  which  raise  the  water  level  encroach  upon  the  freeboard;  and  (3)  surges  may  impose 
sudden  large  loads  upon  miter  gate  operating  machines.  It  should  be  emphasized  that  in  most 
canals,  there  is  no  serious  difficulty  with  surges.  Unless  traffic  density  is  great,  an  existing  surge 
problem  can  quite  easily  be  alleviated.  Nevertheless,  some  consideration  of  the  surge  problem  is 
warranted  in  design,  particularly  for  long,  excavated  canals  controlled  by  high  head  locks. 

б.  Limitations  of  Theoretical  Analysis.  The  theoretical  laws  of  unsteady  flow  in  open  channels 
are  well  established  but,  except  in  a  few  fields  such  as  flood  wave  transmission  and  tidal  hydraulics, 
the  techniques  of  applying  the  laws  to  practical  problems  are  not  yet  well  developed.  Certain 
simplifying  assumptions  which  are  permissible  in  the  study  of  canal  surges  are  entirely  inappro¬ 
priate  for  the  analysis  of  unsteady  flow  phenomena  which  may  occur  in  channels  used  for  other 
purposes.  For  example,  surges  of  large  amplitude  cannot  be  tolerated  in  practical  navigation 
canals  nor  can  the  occurrence  of  high  velocities.  Therefore,  methods  applicable  exclusively  to 
surges  of  small  amplitude  may  be  used.  In  navigation  canals,  the  velocities  of  flow  are  so  low  that 
they  may,  for  some  purposes,  be  neglected.  Limitations  exist  as  to  the  accuracy  which  can  be 
realized,  regardless  of  the  computation  procedures  adopted.  This  is  true  for  several  reasons; 
minor  changes  in  canal  cross  section  produce  partial  reflections  which  cannot  be  evaluated;  traffic 
in  the  canal  has  a  similar  effect;  interference  between  newly  generated  and  reflected  surges  occurs. 
The  latter  problem  would  not  be  a  serious  obstacle  if  the  strength  of  reflected  surges  were  known, 
because  the  effects  of  two  or  more  waves  can  readily  be  combined.  However,  basic  data  as  fo  the 
rate  of  energy  dissipation  in  unsteady  flow  are  unavailable.  For  these  reasons,  and  others  dis¬ 
cussed  in  detail  in  several  of  the  references  listed  at  the  end  of  this  text,  the  results  of  any  analysis 
of  canal  surges  should  be  treated  as  an  index,  comparable  with  similarly  determined  indices  for 
canals  in  successful  operation.  To  avoid  presentation  of  ideas  which  cannot  be  supported  by 
bbaerved  data,  this  chapter  is  confined  to  the  case  of  level  canals  which  are  reasonably  regular  in 
form.  Surges  in  canalized  rivers  are  less  important  than  those  in  canals  because  the  river  widths 
are  generally  much  greater.  This  is  fortunate  because  the  analysis  of  the  effects  of  major  channel 
irregularities  is  very  speculative.  Traffic-induced  surges  in  navigation  canals  are  significant 
mainly  from  the  standpoint  of  vessel  power  requirements  and  are  not  considered  in  this  chapter. 
Surges  which  occur  in  lock  chambers  are  likewise  excluded  as  a  special  problem  in  lock  design  and 
are  discussed  in  chapter  4  of  this  part. 

6-02.  TERMINOLOGY 

Since  the  text  involves  the  repeated  use  of  terms  which  are  not  in  common  use,  it  appears 
advisable  to  introduce  a  few  basic  definitions  at  this  point. 

Surge  is  generally  used  to  describe  the  monoclinical  gravity  wave  of  either  elevation  or  de¬ 
pression  which  is  generated  by  a  sudden  or  rapid  change  in  the  discharge  of  an  open  channel.  A 
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surge  profile  is  governed  by  the  rate  of  change  of  discharge.  In  equalizing  a  navigation  lock  with 
an'adjacent  canal  level,  the  initial  bead  is  great  with  the  maximum  discharge  rate  occurring  as 
soon  as  the  valves  are  fully  opened.  However,  the  head  falls  as  the  equalization  proceeds  and  the 
discharge  gradually  decreases.  Thus,  the  equalization  produces  not  just  one  monoclinical  wave 
but  two,  \vlu<*h  arc  of  opposite  sign;  one  with  a  rather  steep  faci*  whicli  corresponds  to  the  rapid 
initiation  of  discharge,  and  another  with  a  very  gentle  slope,  corresponding  to  the  gradual  reduction 
of  discharge.  The  steep  front  causes  a  departure  from  the  initial  slate  of  rest  in  the  channel  and 
the  gentle  front  of  opposite  sign  causes  a  restoration  to  the  initial  level  and  also  to  the  original 
state  of  Z('ro  How.  (jeiuTaliy  speaking,  somewhat  contrary  to  the  usual  meaning  of  the  term,  tlie 
entire  solitary  wu^•e  of  either  elevation  or  dej)ression  which  is  formed  in  a  canal  as  a  result  of  a  lock 
operation  may  l>e  called  a  surge.  However,  this  solitary  wave  is  really  not  one  surge  but  two  of 
opposite  sign. 

Positive  surrjes  elevate  the  water  surface. 

N'egative  suvjes  depress  the  water  level. 

Celerity  is  the  velocity  of  a  wave  or  clement  thereof,  relative  to  the  velocity  of  flow.  Whoix  a 
wave  is  propagated  through  still  water,  its  absolute  velocity  is  equal  to  the  celerity  but  in  moving 
water  the  absolute  velocity  may  be  cilher  greater  or  less  than  llie  celerity. 

Amplitude  is  the  displacement  of  the  water  level  due  to  a  surge,  measured  vertically  from  the 
original  water  level. 

e-€S.  GENERATION  OF  SURGES 

a.  Smalt  Surges  of  Constant  Amplitude.*’^**  ®'  *  Although  the  profiles  of  the  surges  generated 
by  lock  operations  are  of  varying  amplitude,  it  is  instructive  to  consider  first  the  generation  of 
small  surges  of  constant  amplitude.  Plate  No.  1  presents  the  derivation  of  all  formulas  necessary 
in  the  analysis  of  small  surges  of  constant  amplitude.  With  the  formulas  presented,  it  is  possible 
to  determine  the  celerity,  amplitude,  and  absolute  velocity  of  all  surges,  whether  positive  or  negative, 
propagated  either  upstream  or  downstream.  The  four  possible  cases  are  defined  and  illustrated 
on  Plate  No.  1.  Wlien  the  water  is  initially  at  rest,  the  appropriate  case  may  be  selected  by  making 
K  equal  to  z<  ro,  in  which  case  tlie  celerity  and  (he  absolute  velocity  of  the  surge  are  identical. 
Three  equations  for  the  celerity  are  presented;  the  method  wliich  involves  no  avoidable  approxima¬ 
tions  and,  therefore,  called  the  exact  equation,  an  approximate  equation  named  after  St.  Venant, 
and  another  approximate  equation  named  after  Lagrange.  The  exact  equation  should  be  used  in 
theoretical  studies  and  for  surges  of  great  amplitude.  St.  Venant’s  equation  is  sufficiently  accurate 
for  surges  of  moderate  amplitude.  Lagiange's  equation  may  be  used  for  small  surges.  In  a 
practical  case,  the  three  equations  can  easily  be  tested  by  trial  computation. 

Notation:  The  symbols  used  throughout  the  text  and  on  the  plates  are  listed  at  the  end  of  the 
chapter. 

Sign  Conventions:  In  the  derivations,  shown  on  Plate  No.  1,  all  depths,  velocities,  and  dis¬ 
charges  are  considered  positive;  all  differences  have  been  written  by  deducting  the  lesser  quantity 
from  the  greater  so  the  differences  also  remain  positive.  Careful  study  of  the  definition  sketches 
will  prevent  error  in  sign. 

Continuity  Law:  The  continuity  law  is  used  in  the  form  applicable  to  unsteady  flow  problems. 
The  shaded  area  of  length  S  and  y  represents  the  change  in  cubic  content  of  the  channel  in  one 

^Numbers  shown  thus  refer  to  references  listed  at  end  of  this  chapter.  Wave  mechanics  is  a  problem  in  hydro¬ 
dynamics  which  is  not  w'cll  adapted  to  elementary  treatment.  The  text  has  been  limited  to  readily  comprehensible 
material.  The  references  givcJi  compensate  to  some  extent  for  the  loss  of  exactitude  incurred  in  simplifying  the 
text  to  the  greatest  possible  extent. 


2 


PABT  cxvi,  chapter  6 

March  1949 

sorond.  It  is  evident  that  the  change  in  cubic  content  of  the  channel  whieh  occurs  in  one  second 
must  cqunl  ihe  (lifrerence  hct\vocn  tlie  discharges  into  and  out  of  the  channel. 

Momiittum  Law:  In  applying  the  momentum  law,  a  period  of  one  second  is  considered.  The 
ai)plical)le  force  is  the  ditTerence  in  hydrostatic  pressures  upstream  and  downstream  from  the  surge 
front.  Tiie  'volor-ity  change  is  merely  the  difference  between  the  velocities  prior  to  passage  and  in 
the  wake  of  tlie  surge  front.  The  mass  to  be  used  in  the  momentum  law  must  be  the  mass  which 
is  subject  to  Ihe  V(‘loeity  cl\ange  in  a  period  of  one  second.  It  miglit  appear  that  the  mass  to  be 
used  in  applying  Xewlon’s  Second  Law  corresponds  to  the  volume  but  this  is  erroneous.  A 
surge  traveling  downstream  overtakes  llic  flow  at  a  rate  (S—Vx)  an<l  one  Iraxcling  upstream 
encounters  tlu*  flow  at  a  rate  (Si-Vi).  Thus,  the  mass  to  be  used  in  the  momehluin  equation  is 
wDi/oiS—Vi)  for  surges  traveling  downstream,  an<l  wDi/g{S+Vi)  for  surges  tra  'cling  upstream. 

Form  of  Sfilutioa:  The  momentum  and  continuity  laws  yield  two  simultaneous  equations 
which  may  be  solved  in  any  desired  manner.  Generally,  it  is  <'onvenient  to  eliminate  V2  and  to 
solve  for  as  done  on  Plate  No.  1.  However,  other  solutions  may  be  more  useful  in  some  cases. 
The  expression  for  C  in  terms  of  Di  and  Dz  cannot  be  solved  directly,  since  is  unknown;  suc¬ 
cessive  nppro.ximations  provide  the  most  convenient  solution.  Owing  to  the  notation  adopted, 
the  expression  for  C  is  identical  in  form  for  Cases  I,  II,  III,  and  IV.  It  should  be  noted  carefully 
that  in  Cases  I  and  III  which  deal  with  the  positive  surge,  exceeds  A,  whereas  in  Cases  II  and 
IV,  dealing  with  the  negative  surge,  A  exceeds  A- 

Equations  for  Height  of  Surge:  For  all  four  cases,  an  expression  is  given  for  y,  in  terms  of  Qu 
Qzy  Vx  and  C.  It  should  be  noted  that  it  is  not  necessary  to  use  the  exact  expression  for  C.  When¬ 
ever  St.  Venant’s  and  Lagrange equations  for  celerity  provide  a  sufficiently  accurate  value  of 
A  that  value  of  (7  may  be  used  in  the  equations  for  surge  height. 

6.  Small  Surges  of  Varying  Amplitude.  The  generation  of  small  monoclinical  w’avelets  by 
small  discharge  has  been  illustrated  by  Plate  No.  1.  Any  gradual  variation  of  discharge  can  be 
represented  with  any  desired  degree  of  precision  as  a  succession  of  small  instantaneous  discharge 
changes.  It  is  evident,  therefore,  that  gradually  varying  discharge  is  no  obstacle  to  the  determina¬ 
tion  of  resulting  surge  profiles.  It  is  only  necessary  to  compute  the  celerity  of  each  of  a  series  of 
wavelets  and  to  take  account  of  the  varying  velocity  of  flow.  This  method  of  determining  a  surge 
profile  may  be  called  the  step  method. ***'^-*  The  step  method  may  be  used  not  only  in  determining 
surge  profiles  in  navigation  canals  but  also  other  surges.  For  a  large  and  fairly  rapid  discharge 
change  it  may  not  be  necessary  to  establish  the  surge  profile.  The  exact  equation  for  celerity  may 
then  be  used  directly  to  find  the  surge  amplitude.  The  stationary  hydraulic  jump  is  merely  a 
stopped  finite  surge,  a  phenomenon  which  can  occur  only  when  the  velocity  of  the  oncoming  flow 
in  Case  III  equals  the  celerity  of  the  stopped  surge.**  *  S.  M.  Woodward  and  C.  J.  Pasey  **  discuss 
the  moving  hydraulic  jump  and  show  that  the  absolute  velocity  of  the  jump  is  governed  by  the 
celerity  and  the  velocity  of  flow.  The  jump  may  move  cither  upstream  or  downstream.  The  cose 
of  the  jump  moving  downstream  is  not  covered  on  Plate  No.  1,  but  is  governed  by  principles  identical 
to  those  used  in  that  illustration.  While  these  problems  are  not  directly  pertinent  to  the  study 
of  surges  in  navigation  canals  they  do  present  some  interesting  aspects  which  clarify  the  over-oil 
problem  and  the  concept  of  celerity. and  its  relationship  to  absolute  velocity.  For  the  gradually 
varying  discharges  which  generate  surges  in  navigation  canals,  either  the  step  method  or  a  quick 
approximate  method,  as  described  below,  should  be  used  in  obtaining  surge  profiles.  In  the  step 
method,  each  of  a  series  of  small  discharge  changes  is  treated  by  the  principles  developed  for  surges 
of  constant  amplitude.  However,  because  the  surges  which  occur  in  practical  navigation  canals 
arc  small,  the  effect  of  neglecting  variations  in  the  celerity  with  changes  in  depth  is  also  rather  small. 
Furthermore,  the  celerity  varies  with  the  square  root  of  the  depth  so  the  percentage  variation  in 
celerity  is  about  one-half  as  great  as  the  percentage  variation  in  depth.  The  Lagrange  equation, 
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which  adopted  herein  for  all  celerity  computations^  gives  the  following  variation  of  celerity  with 


depth: 


Depth 

Celerity 

(/i.) 

(Jt.  per  sec.) 

8 

16.1 

10 

18.0 

12 

19.7 

15 

22.0 

20 

25.4 

25 

28.4 

30 

31.1 

40 

35.9 

50 

40.2 

It  is  clear  that  celcritie^^  are  very  high  in  comparison  with  the  velocities  which  can  be  permitted  to 
occur  in  a  navigation  canal.  However,  the  maximum  velocity  of  flow  which  accompanies  the  pas¬ 
sage  of  a  surge  caused  by  a  lock  operation  need  not  bo  very  high  to  cause  some  trouble  in  navigation. 
The  velocity  has  a  very  sudden  onset  which  a  pilot  cannot  anticipate.  Thus  the  velocity  of  flow 
due  to  a  surge,  despite  its  small  magnitude  and  its  relatively  short  duration,  is  much  more  dangerous 
than  a  comparable  velocity  which  is  constant  or  which  arises  gradually  like  a  tidal  current.  Some¬ 
times,  and  this  is  particularly  true  of  deep  canals,  the  velocity  of  flow  is  so  low  in  comparison  with 
the  celerity  that  a  reasonably  correct  surge  profile  can  be  obtained  by  use  of  the  assumptions  that: 

The  ejffeci  of  depth  variation  upon  the  celerity  of  all  elements  of  the  surge  profile  is  negligible. 
On  the  basis  of  this  assumption,  all  elements  of  the  profile  are  assumed  to  be  transmitted  at  the 
celerity  corresponding  to  the  original  depth,  using  the  Lagrange  equation. 

The  effect  of  flaw  velocity  upon  the  absolute  velocity  of  all  elements  of  the  surge  profile  is 
negligible.  On  the  basis  of  this  assumption,  the  celerity  is  treated  as  the  absolute  velocity  of 
all  elements  of  the  surge. 

These  assumptions  and  the  approximate  mctliod  based  thereon  should  not  be  used  unless  the 
limitations  of  the  method  arc  clearly  appreciated.  However,  the  approximate  method  is  very 
convenient  whenever  iis  use  is  permissible.  The  deeper  the  canal  and  the  lower  the  velocities  gen¬ 
erated,  the  smaller  becomes  the  error  in  use  of  the  approximate  method.  Presentation  of  the 
approximate  method  in  this  chapter  is  justified  in  part  by  the  fact  that  the  variations  in  the  surge 
profile,  subsequent  to  generation,  cannot  be  analyzed  with  precision. 

c.  Illustrative  Example  of  Surge  Generation.  To  aid  in  the  use  of  the  step  method  and  the 
approximate  method,  illustrative  examples  have  been  prepared. 

The  following  basic  data  arc  assumed: 


Lock  Chamber  Dimensiori*: 

Length . 

Width . . 

Lift . 

Canal  Dimensions: 

Mean  Depth . 

Surface  Width . 

Timing  Factors: 

Equalization  Period-. 
Valve  Operation  Period 


600  ft.  between  gates. 
110  ft.  between  walls. 
30  ft. 

14  ft. 

200  ft. 

8  minutes. 

1  minute. 


Since  a  canal  of  relatively  shallow  depth;  a  lock  of  rather  high  lift,  and  short  equalization  and  valve 
operation  periods  have  been  selected  for  the  example,  the  comparison  is  unfavorable  to  the  approxi¬ 
mate  mefbod.  That  is,  the  difference  between  the  surge  profiles  indicated  by  the  comparison  is 
about  as  great  as  will  occur  in  any  practical  problem.  Plate  No.  2  shows  the  principal  curves  for 
the  example.  Figure  a  of  that  plate  shows  the  water  level  in  the  chamber  as  a  function  of  time  for 
the  operation  of  equalization  with  the  lower  pool.  This  curve  was  obtained  by  use  of  differential 
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equations  of  discharge  under  falling  head.  In  a  practical  case,  a  comparable  curve  might  he  ob¬ 
tained  from  mo<lcl  tests.  For  simplicity,  ineriial  effects  have  been  ignored  in  deriving  the  curve, 
therefore,  the  typical  “ovcrtravel”  of  levels  is  absent.  Plate  No.  2,  Figure  b  shows  the  discliorgc 
out  of  the  chamber  as  a  function  of  time.  Plate  No.  2,  Figure  c  shows  the  surge  profile,  in  the  canal 
downstream  from  the  loek  at  the  instant  when  equalization  U  completed.  Amplitudes  obtained 
by  both  the.  approximate  method  and  the  step  method  are.  shown.  The  approximate  method  is 
illustrated  in  Table  1  and  the  step  method  in  Table  2. 

Approximate  Method:  Recalling  that  the  approximate  method  is  based  on  the  assumptions  of 
constant  celerity  and  treatment  of  the  celerity  as  an  absolute  velocity,  the  surge  profile  may  be  fully 
determined  by  the  application  of  Lagrange's  equation  and  the  law  of  continuity  to  a  series  of 
arbitrarily  selected  discharges  sufficient  to  define  the  discharge  hydrograph. 

The  celerity  given  by  Lagrange’s  equation  is 
n  C*.^Dj*=21.2  fl.  per  lec. 

Since,  in  the  approximate,  method,  the  absolute,  velocity  of  ail  elements  of  the  surge  is  constant  and 
equal  to  the  celerity,  the  amplitude  of  any  element  of  the  surge  »  ifivan  by  the  equation  for  surge 
amplitude  (Case  I): 


Since  all  data  in  Plate  No,  1  are  based  on  unit  width,  the  discharge  must  be  divided  by  the  surface 
width  of  the  canal  before  using  the  equations.  In  the  present  example,  the  width  is  200  ft.  and  Qi 
and  Vi  are  zero  so: 

ft  -A 

^  (200X21.2)  4240 

The  station  of  each  amplitude  at  the  end  of  the  equalization  period  is  simply  the  product  of  the 
celerity  and  the  travel  time  (in  seconds),  that  is,  the  time  which  has  elapsed  since  generation  of 
amplitude  in  question.  For  a  negative  surge,  generated  when  the  lock  ia  equalized  with  the  upper 
pool,  the  calculations  are  identical,  in  the  appro.ximate  method,  but  the  surge  amplitude  ia  measured 
below  the  original  level.  The  surge  profile  determined  by  the  approximate  method  is  always 
identical  in  shape  to  the  discharge  hydrograph. 

Step  In  the  step  method,  as  in  tho  approximate  method,  a  number  of  discharges  are 

arbitrarily  chosen  from  the  hydrograph.  However,  the  disdiaiges  should  be  quite  closely  spaced 
on  the  hydrograph  so  that  the  heights  of  tho  individual  wavelets  are  small.  In  the  illustrative 
computations,  the  celerity  and  velocity  of  each  wavelet  are  based  upon  the  preexisting  velocity, 
and  depth,  using  the  Lagrange  equation.  It  would  be  more  accurate  to  use  the  St.  Venant  equa¬ 
tion,  particularly  in  the  portion  of  the  computation  in  which  rather  laige  depth  and  velocity  changes 
occur.  However,  the  St.  Venant  equation  requires  repeated  trial  computations  because  tho  celerity 
is  dependent  upon  the  amplitude  which  is  unknown.  The  error  due  to  use  of  the  Lagrange  equation 
is  tolerable  for  the  study  of  navigstipn  canal  surges,  but  it  should  be  recognized  that  an  approximate 
method  is  being  used.  The  columns  of  Table  2  are  filled  in  on  tho  basis  of  tho  following  instructions: 

Col.  1 — Oeveraiion  Time:  Enter  the  Lime  from  the  inception  of  the  occurronce  of  each  of  the 
arbitrarily  assumed  sudden  discharge  changes. 

Col.  2 — Discharge:  Enter  the  discharges  corresponding  to  the  times  shown  in  col.  1. 

CoL  S — Discharge  Increments:  Enter  the  discharge  changes,  using  positive  signs  for  in¬ 
creased  dischoigc  and  negative  signs  for  decreased  dischaige.  (All  discharge  changes  were  con¬ 
sidered  positive  in  Plate  No.  1,  but  in  tabular  computations  it  is  convenient  to  use  plus  and 
minus  signs  to  denote  respectively  increments  and  decrements  of  discharge.) 

Col.  4 — Initial  Depth:  Enter  the  depth  which  obtains  prior  to  the  discharge  change,  Di. 
Note  that  the  initial  depth,  Di  of  any  step  corresponds  to  the  changed  depth  Dt  of  the  preceding 
step. 

Col.  6— Celerity:  Enter  the  celerity  corresponding  to  D|  using  u.  equation. 

Col.  6— Initial  Velocity:  Enter  the  initial  velocity,  V,,  Corresponding  to  Qj  and  Di.  (All 
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diacharges  Bhoivn  on  Plate  No.  1  are  expressed  in  cubic  feet  per  Bccond  per  foot  of  width 
In  this  example,  the  wndth  is  200  ft.,  so  the  diacharpea  per  unit  wddfch  are  determined  by  inoperv 
tion). 

Col,  7 — Svrge  VHocity:  Enter  the  a1)solutc  surge  i^elncily,  w'hich  is-the  sum  of  C  and 
for  Cases  I  and  11,  tlie  difference  between  Cand  for  Cases  HI  and  IV. 

Col,  8 — Amplitude:  Enter  the  amplitude  of  eae)i  eJement  of  the  surge,  which  is  equal  to  the 
discharge  rhan^  divided  by  S\  in  all  cases.  (In  Plate  No.  I,  the  amplitude,  y,  was  always, 
considered  positive  but,  in  tabular  computations  it  is  convenient  to  use  plus  and  minus  signs 
to  denote,  respectively,  inereirienis  and  decrement?^  of  depth). 

Col.  9 — Changed  VepiK:  Enter  the  depth  subsequent  to  the  depth  change,  Dt,  which  is  thvi 
algebraic  sum  of  £>i  and  y. 

Col.  W — Aeeumulated  Amjditude:  Enter  the  accumulated  amplitude  which  is  the  algcbr'fi.^c 
sum  of  ell  previously  occuring  values  of  y. 

Col.  it — TVtiref  Time:  Enter  the  time  which  remains  until  the  end  of  the  discharge  chang 

Col.  — Stations  Enter  the  product  of  tho  absolute  surge  velodty  and  the  travel  time. 


Tablr  1.* — Determination  of  eurge  profile  by  the  approximate  method 
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in  [Mnl<  No.  2,  Figure  c  it  is  sliown  thnt  the  <{ifT(‘ren(0  in  profiles  obtained  by  the  two  methods  is 
small.  As  previously  mentioned,  (he  (lifference  will  he  smaller  in  a  deeper  canal  or  when  the 
lock  lift  is  smaller  in  eoniiuirison  with  the  rnnnl  depth.  This  is  true  because  the  amplitude-depth 
ratio  and  the  vclocity-relerity  ratio  \Aill  lie  smaller.  Case  I  and  Case  II  in  Plate  No.  1  formulas 
are  npjdicnblo  to  the  surges  generated  downstream  from  the  lock;  Case  III  and  Case  IV  formulas 
should  he  used  for  the  surge  generated  upstream  from  the  lock. 

(>  1)4  REFLECTION  CONDITIONS 

Surg<*s  may  he  reflected  jio.siiively,  negativtdy,  or  parlially.  In  (he  latter  ease,  the  reflection 
may  have  eitliii'  a  positive  or  negative  sign. 

<7.  Positive  Reflection.  When  a  surge  or  any  other  wave  form  encounters  a  barrier,  such  as 
a  dea<l  end,  tlie  wave  reverses  its  direction  and  travels  back  upon  itself  without  change  in  sign.  This 
is  lerino<l  positive  reflection.  Thus,  a  positive  wave  is  reflected  as  a  positive  wave,  a  negative  wave 
as  a  negative  wave.  This  phenomenon  is  known  as  positive  reflection.  The  resultant  profile  and 
tin*  velocity  comlitions  in  the  neighborhood  of  tlie  reflecting  surface  may  be  obtained  by  algebraic 
combination  of  the  ofTects  of  the  reflected  and  the  oncoming  waves.*  At  the  reflecting  surface,  the 
ampliludes  of  oncoming  and  reflected  waves  correspond  to  the  same  point  on  the  profile.  Therefore, 
a  grapli.of  the  stage  at  the  reflecting  surface,  plotted  as  a  function  of  time,  is  proportional  to  a 
stnge-trme  graph  at. a  remote  point  in  the  channel  but  the  ordinates  at  the  reflecting  surface  ore 
twice  as  great. 

b.  Negative  Reflection.  When  a  surge  or  other  wave  form  encounters  a  very  large  increase 
in  the  channel  width,  such  as  a  canal  entrance  to  a  sea  or  large  lake,  the  surge  is  reflected  negatively. 
Thus,  a  positive  wave  is  reflected  as  a  negative  wave  and  vice  versa.  There  is  no  simple  mechanical 
analogue  for  negative  reflection,  making  it  difficult  to  visualize;  hut  H.  Rouse®  offers  a  good  explana¬ 
tion  of  (he  phenomenon  and  E.  I.  Brown*  provides  an  example  of  application  in  a  practical  problem. 
As  in  the  case  of  positive  n  fleclion,  the  resultant  profile  in  the  ncighborhoo<t  of  the  reflecting  surface 
may  ho  obtain(‘(i  by  algebraic  oon)bination  of  the  effects  of  oncoming  and  reflected  profiles.  A 
<  iM  ious  phenomenon  occurs  when  a  wave  form  symmetrical  about  its  midpoint  is  reflected  nega¬ 
tively.  Wlum  (he  midordinate  of  the  profile  is  being  reflected,  (he  resultant  profile  throughout  its 
cni  ir<‘  lengtli  is  momentarily  of  zero  amplitude.  An  observed  imperfeet  example  of  this  phenomenon 
is  shown  in  (he  papers  of  F.  W,  Edwards  and  E.  Soucck.“^  In  complete  negative  reflection,  which 
n»quires  a  very  large  expansion,  no  change  in  the  elevation  of  the  water  surface  occurs  at  the  re¬ 
flecting  section. 

c.  Partial  Reflection.  Complete  positive  reflection  corresponds  to  the  dead  end  and  complete 
mgalive  reflection  corresponds  to  the  very  large  sudden  expansion.  These  cases  are  the  defining 
limits  of  infinite  contraction  and  infinite  expansion  respectively.  It  is  apparent  that  finite  expan¬ 
sions  and  contractions  must  correspond  to  intermediate  reflection  conditions.  It  may  be  further 
noted  (hat  no  reflection  oeoiirs  in  a  channel  of  uniform  cross  section  so  such  a  channel  represents  a 
borderlirn*  case  between  partial  positive  and  partial  negative  reflection.  Cr,  which  may  be  called  a 
reflection  coefficient,  is  defined  ns  (he  ratio  of  the  reflected  amplitude  to  the  original  oncoming  ampli- 
( iiile  of  any  wave  form.  Certain  conclusions  may  at  once  be  drawn  as  to  the  effect  of  width  changes 
uptm  (lie  reflection  coefficient.  Lot  n  represent  the  ratio  of  width  of  the  reflecting  channel  to  the 
width  in  which  the  oncoming  wave  is  propagated.  Thus,  if  b  represents  the  width  of  channel  in 
which  the  surge  approaches  the  reflecting  section,  (nb)  is  the  width  of  the  reflecting  section.  The 
same  cic*  pth  is  nssunu'd  in  both  (liannels.**  The  relationship  between  n  and  Cr  is  showm  below. 

When 

n=1.0 . - . - . Cr'=2ero. 

t»  =  zcro . - . . 

n  is  infinite . . . . - . —  — 1.0. 

n>1.0 .  Cr  is  negative. 

>j<1.0 . . Cris  positive. 
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To  satisfy  the  five  known  conditions  concerning  the  vnrintion  in  the  reflection  cooflficient  with  the 
channel  width  ratio,  the  following  expression  may  be  written  by  inspection: 

c 

This  equation  agrees  with  doiivations  arrived  at  by  E,  1.  Brown®  and  H,  Lanib.“  It  should  be 
noted  that  Browm’s  equation  gives  the  resultant  amplitude,  (hat  is  the  sum  of  oncoming  and 
reflecting  amplitude,  not  the  reflected  amplitude  directly.  The  general  problem  of  reflection  and 
transmission  coefficients  for  section  changes  other  than  simple  wucith  changes  is  beyond  the  scope 
of  this  chapter. 

6-05.  STABILITY  AND  LONGEVITYjlOF  SURGES 

A  freshly  generated  surge  profile  can  he  determined  w^ith  sufficient  accuracy  for  practical  needs. 
The  manner  in  which  a  surge  is  reflected  by  a  width  change  can  be  analyzed  without  difficulty. 
However,  the  apparently  simpler  problem  of  determining  the  changes  which  occur  in  the  profile  of 
a  surge  as  it  traverses  a  canal  is  not  on  a  satisfactory  basis  for  practical  applications.  Since  prob¬ 
lems  in  stability  and  longevity  of  surges  are  too  complex  for  treatment  in  this  chapter,  study  of  tlie 
references  is  suggested.  Boussinesq,  ®  ®  has  investigated  the  stable  solitary,  positive,  gravitational 
wave,  his  equation  involves  not  only  amplitude  ami  depth  but  also  the  surface  curvature.  From 
this  study  and  a  similar  one  by  Rayleigh,®  it  is  possible  to  obtain  the  profile  of  the  positive  solitary 
wave  of  stable  form,  that  is,  one  whose  elements  all  have  llic  same  celerity.  It  has  been  suggested,® 
that  the  equations  of  Boussinesq  and  Rayleigh  be  used  for  predicting,  at  least  qualitatively,  the 
trend  of  profile  variation  for  profiles  which  do  not  conTspond  to  the  stable  form.  This  possibility, 
no  doubt,  has  some  application.  Essentially  stable  profiles  have  been  observed  to  travel  many 
miles  with  inappreciable  change,  yet  the  elementary  expressions  for  celerity,  based  only  upon  depth 
and  amplitude,  require  a  positive  solitary  w'ave  to  develop  an  abrupt  front  face  and  an  ever- 
flattening  back  face.  This  characteristic  of  the  elomenttuy  expressions  for  celerity  was  exhibited 
by  the  step  metliod  calculation  of  a  surge  profile  in  the  illustrative  example  of  surge  generation. 
Instead  of  quickly  developing  an  abnipt  front  as  the  computation  indicates,  it  appears  that  the 
surge  develops  a  gently  rounded  profile  of  fairly  stable  form.‘°  Abrupt  fronts,  which  should  be 
clearly  visible,  linvc  not  been  observed  in  navigation  canals.  Furthermore,  it  has  been  shown  ®* 
that  an  abrupt  breaking  surge  of  the  shock  type  can  occur  only  when  the  surge  amplitude  exceeds 
the  dept  For  the  amplitude-depth  ratios  which  occur  in  navigation  canals,  surges  assume  the 
smooth  u  idular  form.  For  these  reasons,  it  may  be  concluded  that  a  solitary  w'ave  generated  by 
lock  operations  develops  a  fairly  stable,  form,  particularly  in  the  case  of  the  positive  wave.  It  is 
known  that  a  solitary  w^ave  of  depression  is  less  stable  than  a  solitary  w  ave  of  elevation.  As  an 
aid  to  the  judgment  in  estimating  the  stability  of  wave  forms,  the  observations  of  J.  Scott  Russell,® 
although  made  more  than  100  years  ago,  are  of  <  onsiderablc  value.  The  problems  of  wave  stability 
are  too  complex  for  complete  treatment  in  this  chapter,  being  w'ithiii  the  province  of  hydrodynamics,*® 
Even  if  a  stable  form  is  assumed,  losses  of  energy,  however  small,  must  occur  because  w^ater  in 
motion  is  involved.  Ultimately,  the  expenditure  of  energy  must  reduce  the  concentration  of  a 
solitary  wave  about  its  mid-point.  No  satisfactory  theoretical  basis  is  available  for  determining 
energy  losses  and  the  corresponding  profile  changes.  Sonic  of  RussclFs  observations  showed  the 
decrease  in  maximum  amplitude  of  a  solitary  wave  as  a  function  of  distance  traversed.^  ®  Additional 
data  arc  needed  to  place  the  problem  of  energy  loss  in  solitary  waves  on  a  sound  basis.  If  the  rate 
of  energy  loss  wtfp  determinahlo,  the  corresponding  pi*ofile  flinnges  could  he  determined  because 
the  potential  and  kinetic  energies  of  waves  can  bi^  <*valujited.''’ The  questions  of  stability  and 
longevity  of  surges  are.  of  considerable  practical  imporliince  in  the  case  of  navigation  canals.  A 
surge  problem  of  greater  magnitude  will  exist  if  the  surges  generated  tend  to  retain  their  strength 
and  integrity  over  long  periods  of  time  than  if  the  surges  tend  to  become  undular,  break  up,  and 
lose  their  strength.  This  is  true  because  the  synebronisra  of  the  effects  of  newly  generated  and 
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reflected  surges  mny  create  severe  conditions.*®  It  is,  therefore,  unfortunate  that  a  t)ractical  basis 
for  estimating  tlic  longevity  of  surges  cannot  he  suggested  at  this  time.  Observation  of  surge 
conditions,  perhaps  purposefully  created,  in  existing  canals  is  at  present  tlie  best  basis  for  cstimat- 
ing  not  only  the  stability  of  surge  profiles  but  the  longevity  of  surges  as  well.  Only  in  this  manner 
can  the  eflTecIs  of  the  minor  partial  reflections,  which  occur  in  unlined  canals  and  numerous 
other  uncertainties,  he  reliably  evaluated. 

6-06.  SURGE  CONTROL 

If  the  costs  of  surge  control  arc  prohihitive  or  if  control  is  otherwise*  impracticable,  the  only 
solution  is  to  anticipate,  evaluate,  nn<l  moke  the  best  design  possible  for  thw  conditions  discussed 
in  paragraph  6-01  a.  If,  however,  some  posit iv(»  action  is  warranted  to  prevent  or  alleviate  a 
surge  problem,  the  following  specific  methods  should  be  oonsidere<l. 

a.  Reduction  in  Equalization  Rate,  It  has  been  shown  that  surge  amplitvides  are,  in  gtneral, 
proportional  to  discharge  rates.  Therefore,  a  surge  con  be  limited  to  any  desired  degree  hy  increas¬ 
ing  the  lock  equalization  period.  In  a  new  design,  this  would  he  done  hy  providing  culverts  of 
low  capacity;  in  an  existing  canal,  the  culverts  wouhl  bo  throttled  in  any  suitalde  manner.  Delay 
to  traffic  is  the  only  objection  to  the  use  of  this  method.  Since  this  solution  is  tlie  only  one  wdiich 
can  be  adopted  at  very  small  cost,  it  should  always  he  consuloretl  very  carefully. 

b.  Increase  in  Channel  Dimensions.  An  increase  in  channel  depth  increases  the  celerity 
and  reduces  the  surge  amplitude;  an  increase  in  w'idth  decreases  the  surge  amplitiKle  without  change 
in  celerity.  In  cither  case,  the  velocity  of  flow  due  to  a  surge  is  reduced  ns  is  also  the  collision 
hazard.  The  cost  is  usually  an  important  objection  to  the  use  of  generous  channel  dimensions  foj* 
surge  control.  Of  course,  if  large  channel  dimensions  are  needed  for  any  otlier  purpose,  the  surge 
problem  is  benefited  incidentally. 

c.  Surge  Control  Reservoirs.  The  following  discussion  of  a  surge  control  reservoir  covers 
only  the  control  of  surges  in  an  upper  pool;  it  is  evident  that  a  similar  up])lication  is  ]>ossil)Ie  for 
equalization  of  the  lower  pool.  The  surge  control  reseivoir  has  small  conduits  or  otlnu*  r(*slric((‘d 
inlets  joining  the  reservoir  and  the  upper  \)ool  Tlic  lock  has  culverts  leading  to  hotli  the  reservoir 
and  to  the  upper  pool  diicctly.  In  raising  the  level  in  the  clinmher,  the  w'ater  is  first  drawn  to  the 
greatest  feasible  extent  from  the  reservoir  hut  the  final  equalization  is  accomplished  with  water 
taken  from  the  iii)i)er  pool.  The  relative  amounts  drawm  from  the  two  sources  depend  upon  the 
reservoir  area.  To  prevent  excessive  equalization  periods,  the  inflow  to  the  lock  can  he  swit<  bed 
from  the  reservoir  to  the  upper  pool  before  the  lock  and  reservoir  has  been  full,y  cqunlizt^d.  The 
small  outlets  conne<*ting  the  reservoir  and  upper  pool  should  he  proportioned  or  controlled  so  as 
lo  just  refill  the  reservoir  prior  to  the  next  lo<  knge.  Thus,  the  rate  of  draft  of  riie  restuwoir  upon 
the  upper  pool  is  very  low  and  the  siirges  due  to  refilling  of  the.  reservoir  arc  negligible.  The 
[>rinripal  objections  to  the  use  of  a  surge  control  re.servoir  are  the  costs  of  the  reservoir  itself  and 
(he  cost  of  the  additional  culverts,  valves,  and  controls.  However,  in  an  important  canal,  the  use 
of  a  surge  control  reservoir  may  he  feasible,  particularly  if  favorable  topographic  conditions  exist. 
A  surge  control  reservoir  was  considered  in  1940  for  the  Panama  Canal  Thiixl  Locks  project.  The 
model  tests  indicated  that  the  reservoir  would  have  been  very  effective. 

d.  Channel  Expansions.  The  phenomenon  of  partial  negative  reflection  at  a  channel  expan¬ 
sion  suggests  the  possibility  that  such  e.vpansions  could  bo  used  to  limit  surge  nmplitiul(‘s.  When 
the  surg<'  encounters  the  expansion,  w'hich  should  he  located  very  luuir  the  lock,  partial  negative 
reflections  will  occur.  Since  the  sign  of  (he  negalividy  r<*flec((‘d  wave  is  opp<»sile  to  that  of  tlu^ 
oncoming  w'avc*,  some  cancellation  w’ill  occur.  It  is  ne<*cssurv  to  coiisich'i*  the  transmitted  (’oni- 
ponent  which  enters  the  e.xpansion  and  its  positive  partial  reflection  from  the  contraction  following 
the  expansion  in  a  complete  anabasis.  The  effect  of  a  channel  expansion  can  be  <lelermined  only 
by  carefully  ooinbining  the  effects  of  all  components  in  accordance  with  the  tlieory  of  partial 
reflections.  It  is  evident,  however,  that  some  beneficial  effects  can  be  obtained  in  this  manner. 
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In  a  favorable  topographic  situation  the  use  of  channel  expansions  might  be  a  feasible  method  of 
surge  control,  although  no  purposeful  application  of  this  method  is  known. 

6-07.  MODEL  TESTS 

Since  the  generation,  propagation,  and  reflection  of  surges  are  essentially  dynamic  and  gravi¬ 
tational  phenomena,  model  tests  based  on  Froude^s  Law  can  be  used  in  the  study  of  surges.  Until 
the  laws  of  energy  loss  in  waves  are  better  known,  some  question  will  exist  as  to  the  accuracy  with 
which  the  longevity  of  surges  is  represented  by  a  model  test.  However,  the  results  are  probably 
as  good  or  better  than  those  obtainable  by  any  other  method.  Distortion  is  permissible  for  some 
purposes,  for  example  studies  of  the  effects  of  channel  expansions.  For  tests  which  must  represent 
the  behavior  of  model  vessels  at  junctions  and  bends,  the  model  should  be  undistorted.  Distortion 
is  also  inadvisable  in  any  study  in  which  the  longevity  of  surges  is  an  important  factor.  Observa¬ 
tional  sensitivity  is  the  most  important  criterion  governing  the  model  size.  The  wave  forms  gen¬ 
erated  by  lock  operations  are  so  long  that  surface  tension  probably  has  little  effect,  even  if  the  model 
is  quite  small.  However,  if  a  model  sui^e  assumes  an  undular  form,  surface  tension  might  affect 
the  celerity.  Stilling  wells  or  other  devices  which  involve  time  lag  cannot  be  used.  Direct  observa- 
tipn  of  water  levels  on  streamlined  staff  g^es  is  feasible ;  good  results  were  obtained  in  a  model  pro¬ 
totype  comparison.’®  The  very  best  sensitivity  which  can  be  expected  from  directly  observed 
gates  when  the  stage  changes  rapidly  is  probably  about  0.002  ft.  Thus,  if  the  protot)rpc  surges 
are  to  be  represented  to  the  nearest  0.1  ft.,  the  scale  ratio,  model  to  prototype,  should  not  be  leas 
than  about  1:50.  In  determining  the  model  scale,  the  characteristics  of  the  phenomena  should  be 
considered  carefully.  Waves  which  have  a  length  of  only  several  inches  in  the  model  are  affected 
by  surface  tension  and  do  not,  therefore,  behave  as  do  corresponding  prototype  waves.  WhUe  the 
number  of  available  model-prototype  comparisons  is  small,  available  data  and  theoretical  con¬ 
siderations  indicate  that  model  studies  of  surge  problems  should  be  quite  reliable  and  trustworthy. 
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